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Abstract: Taking an area of about 2.3X10* km? of southeastern Iran, this study aims to detect and predict 
regional-scale salt-affected lands. Three sets of Landsat images, each set containing 4 images for 1986, 
2000, and 2015 were acquired as the main source of data. Radiometric, atmospheric and cutline blending 
methods were used to improve the quality of images and help better classify salinized land areas under the 
support vector machine method. A set of landscape mettics was also employed to detect the spatial 
pattern of salinized land expansion from 1986 to 2015. Four factors including distance to sea, distance to 
sea water channels, slope, and elevation were identified as the main contributing factors to land 
salinization. These factors were then integrated using the multi-criteria evaluation (MCE) procedure to 
generate land sensitivity map to salinization and also to calibrate the cellular-automata (CA) Markov chain 
(CA-Markov) model for simulation of salt-affected lands up to 2030, 2040 and 2050. The results of this 
study showed a dramatic dispersive expansion of salinized land from 7.7 % to 12.7% of the total study 
area from 1986 to 2015. The majority of ateas prone to salinization and the highest sensitivity of land to 
salinization was found to be in the southeastern parts of the region. The result of the MCE-informed 
CA-Markov model revealed that 20.3% of the study area is likely to be converted to salinized lands by 
2050. The findings of this research provided a view of the magnitude and direction of salinized land 
expansion in a past-to-future time period which should be considered in future land development 
strategies. 
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1 Introduction 


Over the last few decades, land use/land cover (LULC) change has been considerably taken place 
across the entire world and caused substantial adverse impacts on natural resources and 
ecosystem functions (Lambin and Geist, 2008). The impacts of LULC change is heavily 
dependent on the intensity of both natural and anthropogenic driving forces (Meyfroidt et al., 
2013), the affected natural resources and the resilience of ecosystem functions (Wu et al., 2013; 
Xu and Grumbine, 2014; Houghton and Nassikas, 2017). At local scales, the magnitude and 
spatial direction of LULC changes could be properly investigated through field surveying and 
analyzing low altitude aerial photography (Hurskainen and Pellikka, 2004; Benedek and Sziranyi, 
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2009) such as peripheral expansion of an urban area; however, it is a major concern to find out 
how much and to what extent a particular type of LULC conversion occur regionally or globally 
and affect ecological functions (Foley et al., 2005; Lambin and Geist, 2008) such as the loss of 
Arctic sea ice, snow cover, tropical forest lands, continental surface water supplies and the 
expansion of salinized land in arid parts of the world. 

Substantial progress in space-borne remote sensing (RS) techniques through the last five 
decades marked the dawn of a new era for synoptic and periodic observations of the earth's land 
surface at regional levels (Lillesand et al., 2014). Currently, a considerable number of RS image 
capturing systems are capable to collect and transmit huge quantities of data back to earth, 
providing scientists with the opportunities to conduct regional scale mapping (DeFries et al., 
2000; Thenkabail et al., 2009) and to detect changes (Kim et al., 2014; McDowell et al., 2015). A 
large body of these studies took advantage of Landsat imagery data because of systematically 
revisiting and acquiring images from all parts of the earth's surface dating back to 1972, with the 
launch of Landsat 1 (Formerly TIRS-1) and free access to its data archive (Roy et al., 2014; Zhu 
et al., 2015). Furthermore, multi-temporal RS images are considered as the primary data source 
for most LULC change models (Bhatta, 2010). These models could dynamically simulate and 
predict future states of LULC classes in a given region according to the transition rules and the 
associated driving forces which are extracted through analyzing multi-temporal LULC maps 
(Clarke et al., 1997; Eastman, 2012). Under such a context, it is feasible to accurately and 
realistically estimate LULC changes in large spatial and temporal scales, e.g., regional land 
salinization change detection and simulation (Zhou et al., 2012). 

Salinization is ranked as one of the most soil-devastating types of LULC change causing a 
wide spectrum of adverse economic and environmental impacts (Shrivastava and Kumar, 2015). 
Currently, an area of about 1.0x10° hm? (7% of the Earth's land surface) is salinized (Metternicht 
and Zinck, 2003) over the world mainly due to improper management of irrigated agricultural 
practices, low precipitation and high evaporation in arid and semi-arid regions (Rozema and 
Flowers, 2008; Shrivastava and Kumar, 2015). It is estimated that more than 50% of land 
resources and arable regions will be significantly affected by salt by 2050 with an annual 
soil-salinity growth rate of 10% (Jamil et al., 2011). From an economic point of view, soil 
salinization decreases the productivity of land for crop production and threatens social economy 
(Shahbaz and Ashraf, 2013; Shrivastava and Kumar, 2015). Moreover, it could substantially 
influence the global ecological functions and services such as triggering vast sand and dust storms 
in arid and semi-arid regions (Moradi, 2016). Fortunately, RS techniques have facilitated a better 
understanding of the intensity, impacts and the extent of regional scale land salinization 
(Metternicht and Zinck, 2003), especially in sub-Saharan Africa and the Middle East, including 
the study area of this research, where frequent sand and dust storms originated from salinized soil 
lands have become a major problem for people living in this regions and beyond (Moradi, 2016). 
Wu et al. (2014) investigated the expansion of salt land areas caused by improper agricultural 
activities in central Iraq and developed a salinized soil mapping methodology to evaluate 
salinized land expansion from multi-temporal Landsat images. Shalaby and Tateishi (2007) 
applied multi-temporal Landsat images from 1987 to 2001 to detect large-scale salt land 
expansion in the coast of northwestern Egypt, which showed a continual expansion of salt land in 
transitional areas between land and sea water. Lin et al. (2006) and Zhou et al. (2012) attempted 
to simulate the spatial expansion of salt-affected landscapes using LULC change models. These 
studies substantially explored and modeled the effects of driving forces in the expansion of 
salinized land areas. 

In southern Iran, salinized land expansion adversely impacted economy and environment and 
disturbed human life (Moradi, 2016). However, no attempt has yet been made to estimate the 
expansion of large-scale salinized land areas in the area. A national report conducted by Moradi 
(2016) indicated that the southeast part of Iran has become a major source of dust storm in the 
shores of Oman Sea and even Afghanistan. Hence, it is important to map and detect large-scale 
salinized land expansion and to provide a comprehensive picture of salinized land expansion. 
Multi-temporal Landsat imagery data of 1986, 2000, and 2015 and the cellular-automata (CA) 
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Markov chain (CA-Markov) model were used to predict the ongoing future pattern and expansion 
of salinized land areas up to 2030, 2040 and 2050. 


2 Materials and methods 


2.1 Study area 


The study area is approximately 2.3x10* km? located in southeastern Iran (25°03'-26°13'N, 
47°40'-61°58’E; Fig. 1). The altitude of the area is from 0 to 1300 m in the south at the northern 
coast of Oman Sea up to the north. The area is dominated by Indian monsoon in the east and 
Mediterranean climate in the west. The annual mean temperature and mean annual precipitation 
vary significantly across the region. 

Bare lands and rocky surface mountains are the main landscapes in the northern study area and 
salt lands expand across the south and southeast. The coastline of the study area is considered as 
one of the major mangrove habitats of western Asia with an area of about 1500 hm” largely 
covered by Avicennia marina sp. (Fig. 1), supporting a great diversity of marine and terrestrial 
species (Zahed et al., 2010). Although the region is rich in natural resources, temporary 
agriculture, small-scale oyster/shrimp farming activities as well as a small number of residential 
areas, which are sporadically distributed near the coastline, are the main land use types of the 
study area, indicating that this region is poorly planned and managed according to national land 
use planning policies (Moradi, 2016). 
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Fig. 1 Layout of the study area 
2.2 Data 


Landsat imagery was used as the primary source of data to delineate salt land expansion in a 30-a 
time profile. It should be noted that the salinization process in this research was referred to any 
kind of processes which cause the accumulation of soluble salt in surface soils (Hanin et al., 
2016). Such a definition includes visual and spectral changes which were detectable by satellite 
sensors and the users. Due to availability of ancillary data such as GPS-field collected reference 
points and historical LULC maps as well as the magnitude of salt land expansion in the study 
area, Landsat imagery data was selected to explore salinized lands for the years 1986 (TM 
sensor), 2000 (ETM+ sensor), and 2015 (OLI sensor). Moreover, in each year, four Landsat 
images with path and row of 158-42, 157-42, 156-42, and 156-43 were acquired to cover the 
whole extent of the study area (12 tiles in total). Table 1 provided the main characteristics of the 
Landsat images processed in this research. Furthermore, a 30-m spacing digital topographic 
model (DEM) map and detailed digital maps of coastline and sea water channel networks from 
the national project of Makran region development (Moradi, 2016) were employed to produce the 
sensitivity of land to salinization and also to calibrate the LULC change model. 
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Table 1 Research data used in this study 


Sensor Path-row Date Cloud coverage (%) Data type Quality 
156-43 10 June 1986 0 LIT 9 
T 157-42 16 May 1986 0 LIT 9 
156-42 10 June 1986 0 LIT 9 
158-42 8 June 1986 0 LIT o) 
156-43 6 June 2000 0 LIT 9 
re 157-42 7 June 2000 0 LIT o) 
156-42 6 June 2000 0 LIT 9 
158-42 14 June 2000 0 LIT 7 
156-43 26 June 2015 0 LIT 9 
157-42 17 June 2015 0 LIT 9 
oe 156-42 26 June 2015 0 LIT 9 
158-42 8 June 2015 0 LIT 9 


2.3 Image processing and changes detection 


The study area lies within a mosaic of 4 Landsat tiles with different visiting time (different 
number of rows) using different Landsat sensors and pre-processing procedures were conducted 
to prepare the images for classification of salinized land areas. The atmospheric and radiometric 
corrections were simultaneously performed using the FLAASH (fast line-of-sight atmospheric 
analysis of hypercube) commercial module implemented in ENVI v5.3® software (Module, 
2009). FLAASH is a moderate resolution atmospheric transmission radiative transfer-based 
atmospheric correction tool for visible to 3-um spectral range which estimated aerosol/haze-free 
reflectance values based on the spectral behavior of dark pixels observed by Kaufman et al. 
(1997) over water and vegetation covers or employing an automated linear regression method 
over dry soil regions (Cooley et al., 2002; Module, 2009). Input Landsat spectral bands were 
transformed to radiance values in band-interleaved-by-line (BIL) format (Module, 2009) to 
execute FLAASH module. Mid-latitude summer atmospheric and maritime with Kaufman-Tanre 
aerosol retrieval models were also employed to effectively retrieve column water vapor and 
aerosol amount of dark-land pixels over the image scenes (Module, 2009). Finally, using 
additional required parameters from Landsat header files such as flight date and time as well as 
the selected atmospheric correction models, all Landsat tiles were transformed to a common 
atmospheric and radiometric datum. 

In case of image mosaicking, cutline blending (feathering) technique (Herman et al., 2000) 
with edge-to-edge polylines crossing the images between salinized land polygons was utilized to 
smooth the overlapping area of the images. For a detailed description about image feathering 
techniques, please refer to Tutorial (2010). The reflectance values of the lookup table derived 
from the overlapping areas were not significantly different in the images with the equal row 
number of 42 and their data ranges were not balanced through histogram matching process. The 
156-43 images, which contained a large number of dark pixels located over the Oman Sea, 
revealed different spectral values than the above-located 156-42 images in the overlapping 
regions, the spectral values of the 156-43 images were effectively adjusted to the 156-42 image 
scenes to overcome this problem by assigning a linear histogram stretch value of 4%. 

With a set of mosaicked and atmospherically corrected image tiles covering the whole extent of 
the study area, the support vector machine (SVM) technique was employed to identify and map 
salinized land areas for 1986, 2000 and 2015. SVM is a common and well-documented 
supervised classification technique which could obtain accurate classification results especially 
under a limited number of training sample units (Mountrakis et al., 2011). Such ability is of great 
help in LULC classification for large areas where it was not practically possible to collect a broad 
number of sample units. Approximately 1000 pixels related to salinized land areas and 5000 
pixels for the other LULC types were collected through interpretation of a set of False Color 
Composite (FCC) images to train the SVM algorithm. The spectral profile of every training 
sample pixel was intuitively compared with the spectral library of salinized soils collected by 
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Moradi (2016) through extensive field in-situ measurement in the southeastern Iran. 

Spectral signatures of salinized soils against Landsat 8 spectral bands are presented in Figure 2 
providing a better view of training sample units. The grey part of the graph covered all kinds of 
mildly to strongly salinized soils across the study area. Although a variety of definitions of soil 
salinization have been suggested, in line with other remote sensing studies of salinized land 
identification and classification, this paper would rely on the spectral signatures of salinized soils 
detected through interpretation of the spectral behavior of soils to explore salinized lands. The 
premise of this approach is based on the observation of spectral signatures illustrated in Figure 2, 
however, it mainly refers to any kind of processes which caused to the accumulation of soluble 
salt in surface soils (Hanin et al., 2016). 
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Fig. 2 Spectral signature of salinized soil across the study area 


A pixel-based procedure was used to evaluate the accuracy of the classification process based 
on overall accuracy and Kappa coefficient statistics. Reference sample points were equally 
distributed between marginal and interior pixels of salinized land polygons to avoid biased 
accuracy results, especially the effect of positional error (see Lunetta and Lyon, 2004). A small set 
of landscape metrics including the number of patches (NP), patch density (PD), mean patch size 
(MPS), mean nearest neighbor distance (MNN) and largest patch index (LPI). NP and MPS 
metrics were chosen to investigate the degree of dispersive/organic (edge growth) pattern of salt 
land expansion. Spatial continuity of salt land expansion was also explored through quantification 
of PD and MNN. LPI was used to estimate the proportional growth rate of the largest salt land 
patch located in the southeastern study area. Table 2 provided a brief explanation of the selected 
landscape metrics in this research. 


Table 2 Description of the selected landscape metrics (McGarigal and Marks, 1995) 


Landscape metrics Description 

NP >, ; n, total number of patches; n;, the i™ salt land patch. 

PD (plants/hm?) nl(Ax10°); A, total landscape area. 

MPS (hm°) >, (ai / (nx 10°)); ai, the area of i™ patch. 

MNN (m) > (hi/n) ; hi, the nearest edge-to-edge distance of patch i to the nearest neighbor patch. 
LPI (m) max’, (ai / A) ; aj, the largest salt land patch. 


Note: NP, number of patches; PD, patch density; MPS, mean patch size; MNN, mean nearest neighbor distance; LPI, largest patch 
index. 


2.4 CA-Markov model 
We simulated the ongoing probable picture of salt land expansion according to the immediately 
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preceding spatiotemporal changes of salinized land areas from 1986 to 2015. The MCE 
(multi-criteria evaluation)-guided CA-Markov land change model implemented in Idrisi Selva 
(Eastman, 2012) was used to project salinized land expansion up to 2030, 2040 and 2050. 
CA-Markov model was a combined dynamic approach which could realistically simulate LULC 
changes of any given landscape under both biophysical and anthropogenic driving forces 
(Hyandye and Martz, 2017; Palmate et al., 2017), such as prediction of regional-scale salinized 
land expansion (Lin et al., 2006; Zhou et al., 2012). 

In this model, Markov chain estimates the transition probability of a given state (a LULC class) 
to any other states without considering neighborhood effects and spatial dimensions (Arsanjani et 
al., 2013; Brémaud, 2013). A CA model with a standard 5x5 contiguity filter was integrated 
(Eastman, 2012) to spatially allocate the probable LULC changes in desired areas (pixels). 
CA-Markov model systematically produced a set of suitability maps (also known as conditional 
maps; El-Hallaq and Habboub, 2015) to express the probability of a land parcel (pixel) of a given 
LULC class in changing to the targeted LULC class. By using these suitability maps, CA-Markov 
model often fails to incorporate the effects of driving forces and constraints (Mahiny and Clarke, 
2012). However, it was possible for the users to calibrate the model by introducing their own best 
set of suitability maps. In doing so MCE was used as a well-documented linear combination 
method in the process of land suitability analysis to effectively prioritize every pixel of a given 
region based on a set of weighted (w;) and standardized factors (f;) and constraints (C;) through 
the following equation: 


MCE =} fwi] [C:. (1) 


In this research, four transition probability matrices were built using Markov chain model to 
estimate the magnitude of salinized land expansion for 2015 (by the interpretation of changes 
from 1986 to 2000 (calibration phase)), 2030, 2040 and 2050 (by implementing the 1986 and 
2015 as the earlier and later LULC maps, respectively (prediction phase)). Four factors including 
distance to sea, distance to sea water channels, slope, and elevation were identified as the main 
driving forces imposing salinized land expansion through the whole extent of the study area. The 
elevation and distance to sea express the potential effects of saline groundwater of the Oman Sea, 
which controls the upstream freshwater by extensive dam construction (Moradi, 2016). The above 
two factors substantially represent the effects of groundwater which we were not practically able 
to map in the regional scale. The sea water channel represents an anthropogenic driving force 
providing local communities with the opportunity to conduct agricultural activities with saline sea 
waters. The slope is a general factor expressing the sensitivity of soil affected by large-scale 
salinization. All factor maps were transformed to a scale-independent range of 0 to 255 where 
values close to 255 indicate a higher sensitivity of land to salinization using different types of 
fuzzy membership functions (Table 3). All factors were compared in pairs and the corresponding 
weights were extracted by taking the analytical hierarchy process method (Saaty, 2008; Eastman, 
2012). The weighted and standardized factor maps were linearly combined by using Equation 1 to 
produce the final land suitability map. The resulting map of the MCE procedure was referred to 
land sensitivity map to salinization. 


Table 3 Factor maps, fuzzification method and weighting score used to evaluate the sensitivity of land to 
salinization 


Membership function Number of control points 


Factor Reference Weight 
Type Shape c d 
Distance to sea Zhou et al. (2012) Linear MD 0 100,000 0.07 
DEM This study Sigmoidal MD 30 70 0.19 
Distance to channel Zhou et al. (2012) J-shaped MD 1000 30,000 0.51 
Slope This study J-shaped MD 4 10 0.23 


Note: Consistency ratio, 5.9%; Number of comparisons, 6; MD, monotonically decreasing; c and d denote the first and second fuzzy 
set breakpoints, respectively. 
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Having land sensitivity map to salinization, transition probability area matrix of the probable 
salinized land change between 2000 and 2015, as well as by employing the predefined 5x5 
contiguity filter and the 2000 map as the initial land cover map, the CA model was executed to 
simulate the 2015 salt land expansion using relative operating characteristic (ROC) metric 
(Pontius and Schneider, 2001). The resultant simulation map was compared with the reference 
Landsat-derived salinized land areas to evaluate the model performance in the simulation of salt 
land expansion. The prediction phase was also performed to predict salt land expansion up to 
2030, 2040 and 2050 based on the corresponding transition probability area matrices derived 
through Markovian analysis and the 2015 map as the initial land cover map. 


3 Results 


The results of image classification and accuracy assessment are given in Figure 3 and Table 4. 
Accordingly, all maps obtained overall accuracy (OA) and Kappa coefficient (KC) of more than 
0.94 and 0.90, respectively, indicating the acceptable accuracy of the image classification 
procedure. Based on the results, salinized land areas expended more than 1.19x10° hm? between 
1986 and 2015. In 1986, salinized land areas occupied approximately 1.79x10° hm? of the region 
(7.7% of the total area). With an annual growth rate of 0.28%, salinized land areas reached 
2.30x10° hm? in 2000 (9.8% of the total area). The study area was dramatically salt-affected 
between 2000 and 2015 where salinized land areas expanded with an annual growth rate of more 
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Fig. 3 Salinized land areas in 1986, 2000 and 2015 


Table 4 Description of data used for image classification and accuracy assessment 


som ne Accurac 
Mapping Land Area Training sample Validation sample AEE o 
5 2: 
ae type Cn Source No. of points Source No. of points KC OA 
Salinized 1.79 FCC 1000 eae as 700 
1986 FCC E l 0.90 0.94 
Other 21.54 FCC 5000 ye: eran 5000 
interpretation 
Salinized 2.30 FCC 1000 a 1000 
2000 P R 0.93 0,94 
Other 21.03 FCC 5000 S : 5000 
interpretation 
Salinized 2.98 FCC field 1000 Field-collected 2000 
2015 i survey GPS points 0.95 0.98 
FCC field Field-collected j ` 
Other 20.35 RGS 5000 GPS points 8000 


Note: KC, Kappa coefficient; OA, overall accuracy; FCC, false color composite. 
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than 0.29% and reached to 2.98x10° hm? in 2015 (12.7% of the total area). 

Landscape metric quantification provided detailed information about salinized land expansion 
in the study area (Table 5). NP, PD, MPS, and LPI increased and MNN decreased from 1986 to 
2015. The increase in NP indicated the dispersive pattern of salinized land expansion from 1986 
to 2015 in which new salt-affected land patches were developed and expanded. Under such a 
pattern, the pre-existed and newly developed salinized land patches grew rapidly resulted in a 
substantial increase in MPS. The decreasing MNN values confirmed the dispersive pattern of 
salinized land expansion. As represented in Figure 3, the dispersive pattern of salt land expansion 
was limited to specific regions and the PD values increased from 1.01 to 1.12. The largest 
salinized land patch, located in southeastern study area, also experienced a considerable area 
expansion and merged to large sized patches and the LPI value increased dramatically from 1.35 
to 11.85 (777.0% increase) from 1986 to 2015. The results of landscape matric quantification are 
given in Table 5. 


Table 5 Temporal landscape metric values of salinized land expansion in 1986, 2000 and 2015 


Mapping year Change 

Landscape metric 

1986 2000 2015 Increase % 
NP 23,570 24,231 25,386 1816.00 Wall 
PD (plants/m?) 1.01 1.05 1.12 0.11 10.8 
MPS (hm?) oS 13.70 17.80 8.47 90.8 
MNN (m) 128 125 121 -7.00 -5.4 
LPI (%) 1.35 3.11 11.85 10.50 777.0 


The MCE procedure was used to quantitatively assess the sensitivity of land to explore the 
salinization probability and model calibration. Factor maps of land sensitivity analysis including 
distance to sea, distance to sea water channels, slope and elevation as well as the corresponding 
standardized maps were illustrated in Figure 4. The map of land sensitivity to salinization was 
finally produced through a linear combination of the weights and the standardized maps. As 
illustrated in Figure 5. The sensitivity of land to salinization varied significantly from the south to 
north. The highest sensitivity scores were detected in the south, southwest, and southeast of the 
study area where salt land areas had been historically expanded. Conversely, sensitivity scores 
decreased with the increase in distance to sea, distance to sea water channels, and elevation (up to 
1300 m a.s.1.) towards the north part of the study area. 

The result of model calibration was graphically illustrated in Figure 6. Black colored regions 
with an area of approximately 2.6x10° hm? (88% of total salinized land area in 2015) denoted that 
the salt-affected regions were accurately predicted. The error areas, which were mistakenly 
classified through the modeling procedure, were categorized into two classes including salinized 
land mapped, but not predicted (area in green, 1.64x10* hm?; 4.6% of the predicted salinized land 
areas; Fig. 5) and salinized land predicted, but not mapped (area in red, 1.21x10*+ hm?; 3.4% of 
the predicted salinized land areas; Fig. 5). The former class was mainly distributed in the eastern 
study area while the latter was located in central and western parts. The salinized land areas were 
estimated being 2.0x10° hm?, approximately 7330 hm? less than the actual salinized land 
expansion of 2015, by the 1986-2000 transition area matrix generated through the calibration 
process. The validity of the MCE-informed CA-Markov model was assessed by comparing the 
projected and actual salt-affected land generated for 2015. The ROC for the salinized land area 
was 0.74 indicating a good agreement between observed and actual salinized land areas. 

The CA-Markov prediction phase was conducted to simulate the salinized land expansion up to 
2030, 2040 and 2050 (Fig. 7). The salt-affected land areas for 2030, 2040 and 2050 were 
estimated to be around 3.7x10°, 4.1x10° and 4.7x10° hm? (15.7%, 17.8%, and 20.3% of the total 
area), respectively. In total, it was estimated that the study area has been and will face more than 
3.0x10° hm? (12.6% of the total area) of salinization from 1986 to 2050. 
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Fig. 4 Evaluation factors (left column) and the resulting fuzzy standardized maps (right column) 
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Fig.5 Salinization land sensitivity derived from multi-criteria evaluation procedure 
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Fig. 6 Error salinized land map of the 2015 model simulation 
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Fig. 7 Growth trend of salinized land areas in 2015 and up to 2030, 2040 and 2050 


4 Discussion 


In southeastern Iran, most ecosystem functions and processes had crossed their ecological 
thresholds and were deeply threatening human life and health in a variety of ways, such as the 
severe water scarcity (Foltz, 2002). Hence, land use managers and decision makers have been 
evaluating the suitability of natural-land resources to migrate settlers from environmentally 
devastated ecosystems of central Iran to these regions (Moradi, 2016) especially to the Oman Sea 
shorelines, the study area. However, this area also faced large-scale soil salinization due to the 
direct and indirect anthropogenic driving forces. A high degree of evaporation, extensive 
freshwater diversions and upstream blocking by dams have substantially decreased the humidity 
of the salinized soils (Moradi, 2016). Frequent storms and strong winds carried a considerable 
amount of dried salinized dust passing over the area. According to the USA Today, storms 
originated from this region made Zabol the most polluted city of the world in 2015 (Bacon, 
2016), highlighting the magnitude impact of the salinized land expansion in this region. 

Processed images successfully revealed the magnitude and spatial trend of salt land expansion 
in the southeastern Iran. The results showed that the major sources of soil salinization were in the 
southeast (including the largest salinized land patch) and west parts of the study area. Temporal 
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analysis of landscape metric values also suggested a dispersive pattern of salinized land 
expansion since1986, with the increase in NP and MPS as well as the decrease in MNN values. 
Although following an outward-oriented growth pattern, salt-affected areas were significantly 
limited to the areas near sea shoreline and water channels in low latitude and gentle slope lands 
with increased PD values. The expansion of the largest salinized land patch in the southeastern 
study area was investigated through quantification of LPI metric. A substantial increase in the LPI 
values, particularly during the period 2000-2015, was due to the merge of the two largest 
salinized patches. 

A further attempt was made to quantitatively explore the sensitivity of land to salinization in 
the study area. The assumption behind this effort was based on the identification of those factors 
which control regional scale soil salinization in wide time spans, i.e., over decades. Zhou et al. 
(2012) found that the driving forces imposing salinization, both biophysical and human-induced 
factors, were spatiotemporally dynamic at regional scale. The researches by Lin et al. (2006) and 
Zhou et al. (2012) revealed the ability of CA-Markov model in the simulation and prediction of 
salinized land expansion. This model was employed to project the ongoing salinized land 
expansion for 2030, 2040 and 2050. The results of this research indicate that the southern Iran is 
experiencing an accelerated rate of salt land expansion. According to the past changes delineated 
in this research, salinized lands will be expanded significantly and will likely occupy the majority 
of plain region. This expansion will form a series of environmental crises. For instance, a very 
large salt patch which was predicted to growth in the southeastern study area would have a high 
capability to trigger regional dust storms. Although the results of this research provided a holistic 
view into the magnitude and direction of salt land expansion, further attempts need to increase the 
reliability of these predictions. In case of criteria selection, one should be careful about the 
temporal variability of driving forces. As Zhou et al. (2012) outlined, driving forces contributing 
to salt land expansion are highly dynamic and have variable properties which can change even 
annually. Keeping this matter in mind, it is of high importance to take into account the probability 
of changes in driving forces as well as the growing speed of salinization for the prediction node 
and in turn for adopting proper land use planning strategies. 


5 Conclusions 


This study evaluated and modeled the spatial expansion of salt lands in southeastern Iran. The 
results showed that salt lands have been growing across the entire plain regions of the study area. 
More importantly, it is predicted that such an expansion will lead to the generation of a very large 
salt patch in the southeastern part with probably significant regional destructive effects. The 
analysis carried out in this research were all based on the spatial analysis of salt lands quantified 
from processing Landsat 8 imagery data in three time steps. Furthermore, the predictions were 
supported by a small set of highly dynamic variables. Despite the successful application of the 
simulation model, it is recommended to restrict the area under investigation and shorten the 
prediction spans to better scrutinize how and to what degree the salinized lands are expanding. 
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